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Abstract. 

A parametrization of octupole plus quadrupole deformation, in terms of 
intrinsic variables defined in the rest frame of the overall tensor of inertia, 
is presented and discussed. The model is valid for situations close to the 
axial symmetry, but non axial deformation parameters are not frozen to 
zero. The properties of the octupole excitations in the deformed Thorium 
isotopes ^■^^Th, ^^^Th are interpreted in the frame of this model. A tentative 
interpretation of octupole oscillations in nuclei close to the X(5) symmetry, 
in terms of an exactly separable potential, is also discussed. 



1 Introduction 

We present here a formalism to describe the simultaneous octupole and 
quadrupole deformations of the nuclear surface, close to but not coincident with 
the axial symmetry limit, in the frame of Bohr hydrodynamical model. This 
scheme has been developed in order to discuss phase-transition phenomena in 
the octupole degree of freedom, a subject which - as it will be explained in the 
following Section|2]- appeared as a natural development of the experimental and 
theoretical researches of the group of Florence in the recent years. 

Sections |3]and [^describe the basis of this model and its application to the 
phase transition between octupole oscillations around a permanently deformed 
reflection-symmetric shape and rigid rotation of a reflection-asymmetric rotor. 
These parts summarize the results contained in a recent paper in the Physical 
Review C [1]. Finally, some new results are reported in the Section |5] while 
in Section |6] a simple model, based on a separable potential, is discussed and 
used to tentatively interpret the octupole oscillations in nuclei close to the X(5) 
critical point. 
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2 Some history 

Our interest on the octupole excitations dates at least from the late '90s, when 
we have searched and found evidence for two-octupole-phonon excitations in 
^^^Nd and ^"'^Sm [2] and we have long searched, but not found, three-octupole- 
phonon excitations in ^^^Gd [3]. 

Later, we have been impressed by lachello results [4, 5] on the symmetries 
at the critical point, such as E(5) and X(5) (other symmetries, like Y(5) and 
Z(5), have been proposed later [6, 7]). The first example of X(5) symmetry 
was identified in ^^^Sm [8], a well known transitional nucleus. We had some 
experience of another transitional region, that of Mo and Tc isotopes: so, we 
could identify another possible X(5) candidate [9] in the ^'^''Mo isotope ^ We 
also observed that the X(5) model was able to account for not only the ground- 
state band, but also the excited 7 bands, both in ^°*Mo and in ^^^Sm. This 
fact is not trivial, because the negative parity bands, present at low excitation in 
^^^Sm, show a very different behaviour. We shall come back to this point in the 
following. 

We also noted some regularities in the proton and neutron numbers of these 
two X(5) nuclei. In fact, ^°^Mo has Z = 42, N = 62, ^^"^Sm has Z = 62, N = 
90: one could suspect that another phase transition takes place also for Z w 90. 

Actually, the standard indicators of quadrupole collectivity, and in particu- 
lar the ratio R2 = iJ(4+)/iJ(2+) give indication of a phase transition in the 
Ra {Z = 88) and Th {Z = 90) isotopic chains. Heavier isotopes have rota- 
tional character, while the lighter ones appear to be vibrational (or non collec- 
tive). Moreover, the energies of positive parity levels in the ground-state band 
of ^^^Th and ^^^Ra show an impressive agreement with the predictions of the 
X(5) model [10, 11]. However, in addition to the positive parity band, these 
nuclei possess a odd- J, negative parity band which starts with a 1~ level ly- 
ing slightly above the 2+ and merges with the positive-parity one at J « 5 or 
6. This means that octupole degrees of freedom are important in these nuclei, 
and their effect must be considered when the evolution of the nuclear shape is 
followed along the isotopic chain. We observe that ^^"^Th and heavier isotopes 
of Th give evidence of octupole vibrations combined with the rotation of a de- 
formed (but reflection symmetric) core: rotational-like bands are built over the 
one-octupole-phonon states with different values of K (the angular-momentum 
component along the approximate symmetry axis). In lighter isotopes, the exci- 
tation energy of the 1~ band head of the K'^ = 0~ band decreases well below 
those of the other octupole bands, and approaches the rotational-energy value. 
Eventually, one observes an alternate-parity band which approaches (but never 
reaches) the behaviour of a rigid, rotation asymmetric rotor, typical of binary 
asymmetric molecules. In conclusion, a phase transition for the octupole de- 

' Level energies and 7 branching ratios in ^"''Mo are in excellent agreement with the X(5) model. 
Later measurements [12] of mean lives showed, however, that B(E2) values are at variance with the 
model predictions. 
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grees of freedom seems to take place in the Th and Ra isotopes, not far from the 
critical point of the quadrupole one, but in the opposite direction: with octupole 
vibration where the quadrupole deformation is stable, and vice-versa. We can 
also observe that only the = 0~ band seems to change substantially along 
the isotopic chain. This means that the relevant quadrupole and octupole degrees 
of freedom are those related to axially synmietric, or quasi-axially-synraietric, 
deformations. 

At this point, one needs a theoretical scheme able to describe the evolution 
of quadrupole and octupole deformations, close to the axial symmetry, across 
the critical point of the phase transition. 

In the frame of the algebraic approach, new developments of the Interacting 
Vector Boson Model have been reported at this workshop [13]. The spdf exten- 
sion of the IBM [14] is also able to account for negative-parity excitations. This 
model, however, is not the most suitable to describe the phase transitions. 

In the frame of the geometrical approach, the most complete treatment is the 
one proposed by Donner and Greiner [15]. They use the Bohr intrinsic frame, 
defined by the symmetry axes of the quadrupole, and refer to this frame the seven 
amplitudes of the octupole mode (the overall tensor of inertia being not diagonal 
in this reference frame). This approach, however, is useful only when all the 
octupole amplitudes are small compared to the quadrupole deformation. Other 
models [16-18] frozen to zero part of the dynamical variables, and are usually 
limited to the axially symmetric case: what is probably enough in the presence 
of a stable quadrupole deformation, but not necessarily when the quadrupole 
deformation approaches zero. 

Finally, we have noted a very interesting work by Wexler and Dussel [19], 
which shows that it is possible to define an intrinsic reference frame for the 
octupole mode, in which the tensor of inertia is diagonal. In this frame, the seven 
octupole amplitudes can be parametrized in terms of four intrinsic variables (as 
the five quadrupole amplitudes are parametrized in term of (32 and 72 in the Bohr 
representation). This approach, however, is valid for the octupole mode alone. 

In conclusion: in order to describe quadrupole and octupole deformation, 
close to but not coincident with the axial symmetry and with the octupole ampli- 
tude not necessarily small compared with the quadrupole one, we were obliged 
to introduce our own parametrization scheme. 

3 The basis of the model 

We must start, as usual, from the definition of the nuclear surface in polar coor- 
dinates, in the laboratory frame: 

A=2,3 iJ,= -\,\ 

We limit our model space to quadrupole and octupole, and we assume that the 
deformation is small enough not to need to introduce a monopole or dipole term 
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to keep constant the nuclear volume and the center-of-mass position. For the 
expression of the kinetic energy, we use that of the Bohr hydrodynamical model, 
T = i J2\ B\ |d|j'^^ p . Then, the laboratory ampUtudes a^^' are expressed 
in terms of the corresponding amplitudes in a proper intrinsic frame and of the 
Euler angles 9i, 02, 03 defining the orientation of the intrinsic frame in the lab. 



Bxa 



(A) 



(2) 



where the D^^^ are Wigner matrices. Note that, in order to simplify the no- 
tations, the inertia parameter B\ has been included in our definition of a^^\ 
In terms of the new variables, the expression of the kinetic energy splits in 
three parts: a vibrational term Tvib = 5 X^a /i I^m'^ 1^' ^ rotational term 
Trot = I J2k k' J^kk'QkQk'f quadratic in the intrinsic components qk of the an- 
gular velocity; and a coupling term Tcoup, not present in the Bohr model for 
quadrupole motion: T^oup = [q'^^^ ® [a^^^ (g) a*^)]^^)]^,"^ For the non 

diagonal components of the tensor of inertia, the hydrodynamical model gives 



J13 + i J23 
J23 



(2) 



'^C2(A)Im[o 



W 6?, n W 



X 



(2) 



(3) 



where C2(2) = -V21, C2(3) = \/l26. It is relatively easy to find a 
parametrization which automatically sets to zero the non-diagonal products 
of inertia, as far as quadrupole and octupole are treated separately. For the 



quadrupole, it is the classical one by Bohr, Og 



(2) _ 



,(2) 



/32COS72, al 



(2) 





^1/2 (32 sin 72. For the octupole, we assume a parametrization similar 



to that of Wexler and Dussel [19]: 



4^) = 

a?' = 
a?) = 



cos 73 

-(5/2) (X + iY) sin 73 

(33 sin 73 
X COS73 + (a/T5/2) sin 73 

,(A) 



(4) 



i Y COS73 — (vl5/2) sin 73 



In both cases, the aji amplitudes with = 2 are real, while the = 1 terms 
are either zero or small of the second order, if we consider small of the first order 
the other non-axial amplitudes. 

The problem arises when both quadrupole and octupole terms are present, 
since the principal axes of the quadrupole and of the octupole tensor of inertia 
do not necessarily coincide. 

We want to use, as in the Bohr model, an intrinsic frame defined by the 
principal axes of the tensor of inertia. We must therefore impose the conditions 
Jij = for i ^ j. This is a set of 3 non linear equations, but we can linearize 
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them if we assume that non axial amplitudes are small in comparison with the 
axial ones. One obtains, up to the first order in the small amplitudes, 

Ji2 = -2 V6 (/32lm ^ + Vd Pslm a^^^ ) = 

Ji3 + i^23 = V6 (/32af^ + V2/33af^) =0 (5) 

which are automatically verified if we put 

a'^i'' = -ciV^Pa ijic + <c) = ci/32 {qc + <c) 

(6) 

Im a~2^ = -C2\/5 Im Cj^'' = 02^2 

with the new parameters r/c, Cc and (_c small of the first order The factors ci, C2 
are - at the moment - arbitrary functions of P2 and /Sa . 

Now, we can introduce the matrix of inertia Q, which relates the classical 
kinetic energy to the components Uj of the generaUzed velocity vector: 

7" = ^ ^ GikUiUk (7) 

where, in our case, u = {(32, 72, /^s, 73, X, Y, ?), (, qi, 92, 93}- The 
matrix Q and its determinant G play an important role in the quantization of the 
kinetic energy, according to the Pauli procedure [20] . 

If we want that the results of the Bohr model be obtained at the limit of 
P3 ^ /32, we must verify that this happens, first of all, for the determinant G. 
Now, the simplest possible choice of constant factors ci ad C2 in the Eq.|6]would 
give G cx /Jg^, while in the Bohr model it should be c>c /3|. A better choice 
would be 

ci = iPl + 2/3|)-V2 C2 = {Pi + 5/3|)-V2 (8) 

With this choice, G has the correct limit when (3^ 0. Moreover, all non 
diagonal elements of Q involving either (32 or (3^ and one of the derivative of the 
other intrinsic amplitudes or 53 are exactly zero. Other non diagonal elements 
are small at least of the first order and it is possible to show [1] that they have 
negligible effects, apart from the elements of the last line and column. The latter, 
in fact, are also small of the first order, but must be compared with the diagonal 
element jTs, which is small of the second order 

Relevant terms of the matrix Q are shown in the Table [2 It can be observed 
that the variables 72 and 73 appear only in the combination 7 = \/5 72 — 73. 
In the place of 72 and 73 it will be convenient to use the new variable 7 and the 
orthogonal combination 70 — cq (/3|72 + V5/3373), entering in the expression of 
G. Once more, the factor co is - at the moment - an arbitrary function of (32 and 
(3^. We can observe that 70 is a measure of the triaxiality of the overall tensor of 
inertia. In fact, if 70 — Q, Ji — J2 (apart from terms of the second order) and 
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Table 1 . The matrix of inertia Q : leading terms and relevant first-order terms. Other 
first-order terms are indicated with the symbol ~0. Here 7 — \/572 — 73, J7i = 3(/3| + 
2/3|) -f 2V3(/3|72 + ^/3|73); J2 = 3(/3| + 2/3|) - 2^{(3h'2 + V50h3); and 

J3 = 4(/3|7| + Phi) + 18{X^ + + 2(77^ + C') + 8C'. 
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the tensor of inertia is axially symmetric, even if the nuclear surface is not. In 
some way, 70 plays a role similar to that of 72 in the pure quadrupole case. 

Now, to characterize the different degrees of freedom with respect to the 
angular momentum projection L3 and to the parity, it is convenient to perform 
another change of representation: 



X 


= w sin 1? 


Y 


V 


= vsm(fi 


c 




= u sin X 


7 


70 
Co 


= /o(/32,/33) Mo 





wcos'd 

VCOS(fi 



^/32 + 5/31/ ,92/33) M cos X 



(9) 

Choosing fo{p2, Ps) = \/ Pi + ^P'i one obtains for the determinant of Q 

G = Detg ^ 2304 (/3| + 2/3|)^ ^2 ^2 _ 

At this point, the matrix of inertia takes a very simple form. The relevant non- 
diagonal terms are limited to the smaU sub-matrix involving and the time 
derivatives of (p, x and 1?, which can be diagonalized easily. To this purpose, we 
evaluate the intrinsic angular momentum component L3 = dT/ dqz and also the 
conjugate moments p^, pj^, of the angular variables Xi 
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^ 2u^ (X + 2<Z3) 
= 2w^ {-d + 3(73) 



(11) 



(12) 



One can immediately observe that when uq — (i.e. when J7i = 93 — ^ 00 
unless — = Pip + 2p^ + 3p^. This means that ft is the component of 
the intrinsic angular momentum along the axis 3, which survives also if this is 
an axial-symmetry axis for the overall tensor of inertia. If we assume that the 
potential energy does not depend on (p, x, or it is easy to realize that p^p, p^ 
and ps are quantized and are integer multiples of H: the degrees of freedom 
corresponding to the pairs of variables (w, ip), {u, x), (w, z?) correspond to 1, 
2 and 3 units of angular momentum along the intrinsic axis 3. This is obvious in 
the latter case, which concerns the octupole variables aj_3, but not for the other 
two, in which quadrupole and octupole variables are mixed together It is also 
possible to show [1] that these variables carry negative parity. 

It remains to consider the variable uq. If we assume that the differential 
equation for uq can be decoupled from the others, this equation takes the form 




uo 



d 

dun 



U{uo)] 




(t>iuo) - (13) 



where we have put Kq = L3 — ft. The variable uo can take positive as well 
as negative values. The condition of continuity for the wavefunction 0(7Xo) and 
its derivative at uo = imposes that Kq — 2nu(,, with n^g integer. We can 
conclude that the degree of freedom associated to the variable uq carries two 
units of angular momentum along the intrinsic axis 3, and it is possible to show 
that it carries positive parity. The inverse of the matrix Q turns out to be diagonal 
(at the relevant order) in the space of momenta conjugate to the variables defined 
in the Eq. |9]and of the angular momentum components Li, L2 and Kq. The 
Table |2]shows a more general form of this matrix, with the variable uq replaced 
by 70 = uo/R{p2,p3), where i? = l/[co(/32, /33)/o(/32, /^s)] (seeEq.©. A 
more formal derivation of these results, involving the derivatives of the Euler 
angles, can be found in the Appendix C of ref. [1]. 



4 The axial octupole mode with stable quadrupole deformation 

This is the simplest case in which the properties of octupole excitation can be 
followed from the limit of harmonic oscillations around the reflection symmetric 
core to the opposite hmit of stable octupole deformation. A detailed discussion 
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Table 2. The matrix Q ^ in the space of momenta (only the leading terms are shown). 



Here, Ko — — n, with Q — p^p + 2p^ + 3p^, and 70 = ito/ R. 
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of this subject can be found in the ref. [1], and a few prehminary resuhs have 
been reported at two earUer Conferences [10, 11]. Here, we only summarize 
these results. 

A preliminary comment is in order. The properties of the quadrupole vi- 
brations around an axially deformed core are better described [23] with respect 

.... (2) r (2) (2) _(2) , . , 1 

to the intrinsic parameters , oog = Oq — ag than in terms of the Bohr 

parameters j32 and 72. We have seen that the parameter uq defined in the eq.|9] 

(2) 

plays, in our treatment, a role analogous to that of in the pure quadrupole 
Hamiltonian. It appears reasonable, therefore, to use it as a dynamical variable 
instead of defining an angle variable similar to the 72 of the quadrupole case. 
Therefore, we can use the form of the matrix given in Table|2] with R — 1, 
to derive the differential equation for fi^, according to the Pauli prescriptions (in 
doing this, we assume decoupling of the motion from the small-amplitude 
oscillations in all other degrees of freedom). One obtains 



d2-0(x) 2x d'0(x) 
dx^ 1 + Ax 



6(1 + a;2) 



v{x) 



V'(x) = (14) 



where x — ^/2(3z/P2, while e are the potential energy and the energy 

eigenvalue in a proper energy unit, and ■)/'(— 2:) = {—)'^ %l}{x). As for the potential 
v{x), we have considered two simple cases: a quadratic expression v = \cx^ 
or a critical (square-well) potential, as in the X(5) model: v{x) — for \x\ < b 
and ~ +00 for |a;| > b. In both cases, the model has one free parameter (c 
or b) to be adjusted to fit the experimental data. In the figure the energies 
of positive and negative parity levels of ^^^Th and ^^'^Th are compared with 
different model predictions. The former turns out to be close to the results we 
obtain for a critical-point potential, while the latter is closer to those obtained 
with a quadratic potential. 



Description of the Nuclear Octupole and Quadrupole Deformation 9 




5 1015 20 5 1015 20 

J J 



Figure 1. (From ref. [1]). Excitation energies for states of positive parity (circles) and 
negative parity (triangles), in units of E{2'^), for the ground-state band of ^^*Th and 
228jj^ Theoretical curves: a - rigid rotor; b {b') - present model with critical potential, 
fitted on the 1^ state (on the high-spin states); c - present model with harmonic potential. 



5 Latest developments 

The results reported until now are, apart from some limited extensions, those 
contained in our recent paper in Physical Review C [ 1 ] . From now on, we move 
to the region where not only f]^ but also is allowed to vary, in principle down 
to zero: this is still a partially unexplored land, were work is in progress and 
results might be subject to change. Eventually, it will be useful to move from 
the Cartesian representation in terms of (^2 , /Sa to a polar one, in terms of the new 
variables /3, S with /32 — /3cos5, /Js — f3 sin S as in the work by Minkov [21], 
but for the moment we will continue to use the Cartesian representation. 

For consistency with the well established results on the pure quadrupole case, 
we need that our results converge to those of the Bohr model in the limit of small 
octupole deformation. As we are going to see, this requirement will put some 
restrictions on the possible models of octupole plus quadrupole excitations. E.g., 
the form of the matrix Q we have used in the case of permanent quadrupole de- 
formation is not suitable in the present case, as its determinant G, at the limit of 
small Pz, is proportional to while it should be proportional to /3f according 
to the Bohr model. It is easy to recognize that the responsibility for the disagree- 
ment can be attributed to the choice of uq as a dynamical variable. Since in the 

(2) 

Bohr model the variable 72 is used in the place of % it is now necessary to 
replace uq with an adimensional variable 70 = uq/ R. We could choose now 
R = P = \//3| + Pi, to obtain the correct hmit G cx /3| for /^a ^ 0. This 
is somewhat better, but still not enough: the equation for P^ obtained with the 
Pauli quantization rule does not converge to the one of Bohr when P^ 0. 

In fact, if we only take into account the dynamical variables P2 and P^ (as- 
suming that their equations can be approximately separated from that of 70, as 
in lachello X(5) model, and from all other dynamical variables), we obtain 
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+ 



-[E-Vo{p2,f33)]- 



3(/3| + 2/3|) 



vI/(/32,/?3) = (15) 



This equation can be somewhat simphfied with the substitution 

*(/?2,/33)=ff-'/'$(/32,/33) 

with g oc G^/^, to obtain 
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2/31), and one obtains 



t^(/32,/33) = - 



5/?| 



16/3|/3g + 14/3| 
(/3| + /3|)(/3i + /3|)' ■ 02 



^ --2 for/33 







(19) 



A similar calculation is possible for the pure quadrupole case, starting from the 
Bohr expression of G oc /32, but the result is v{(32) = — 2//32. Even if the Umit 
of G for 6s converge to the corresponding one of the Bohr model, this is 
not necessarily true for v. 

However, it is easy to realize that v(/32, /^s) converges to v{f32) of the Bohr 
model if the first and second partial derivatives of g with respect to fS^ tend to 
zero when /33 — » 0. A possible choice of variables leading to this result corre- 
sponds to keeping cq = \^02 + 2/31 (instead of -\//32 + 5/31 ) in the definition 
of 7o. In this case one obtains 



9(02, P3) 



(/3| 



2\2 



01 + 5/332 



/3| 1 + 8 (/33/ /32)" + 



J2 -r t./^3 

and the first and second derivative of g with respect to vanish for ps 



(20) 
0. 



6 Specific model for quadrupole-octupole oscillations 

We now consider the case of simultaneous quadrupole-octupole oscillations, and 
in particular the quadrupole motion corresponding to the critical point of phase 
transition described by the X(5) model. At the hmit of small amplitude for the 
octupole oscillations, we should therefore obtain the same results of X(5). The 
point is that the octupole amplitude must be, at this limit, small compared to 
the quadrupole amplitude, which, in turn, can become zero. It is convenient, 
therefore, to use the new variables /3, 5 defined by 



02= COS 5 



= P sin 5 



(21) 
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and reach the limit by confining 5 to very small values by a proper potential 
term. It is clear that also the amplitudes v, u, w must be small compared to 
/3, as well as 70 compared to 1 . At the moment, however, we will forget their 
presence and only discuss the Schrodinger equation involving the variables f3, S 
and the Euler angles. With this ansatz, and assuming that the potential energy 
has the form {H^ /2) V{P, 5), for = we obtain 




d_ 

dp 



ld_ 
'gdS 



-V- 



/32(l + sin" S) 



^iP,S)=0 (22) 



where g oc G^/^ ^ p5 ^-^^2 J (1 + 5)] joUovuw,e = 2E/tf, 
and A J = J(J + l)/3. 

The Eq.|22lhas a structure very similar to that of the Bohr equation for pure 
quadrupole motion at the limit close to the axial symmetry, with our parameter 
6 in the place of 72. In the case of X(5) symmetry, the Bohr equation has been 
solved [5] by approximate separation of the variables, substituting the factor 
1//3| with a proper average value in the differential equation for 72. It has 
already been noted, however, that this approach does work for the 7 excited 
bands but not for the negative-parity ones [9]. Now we want to explore some 
alternative procedure which could better account for the experimental data. 

It is convenient to exploit the result of Ea.s I16I17I18I to eliminate in the 
Eq|22|the first-derivative terms, with the substitution 
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givmg 




7 1 2sin2(5 (25 + 12sin2 5 + 8sin^(5) 1 
4 ^ " 



(24) 



(l + sin^ S){l + AsiTi^ Sy 



We can separate the potential term in two parts, one of which depends only on 
f3: V{P, S) = Vo{P) + Vi(S, 13) and search a solution in the form ^'o(/3, 5) = 
V'(/3)1>(5,/3): 



1 



dp 



Aj + l + A' 



/32 



(25) 



1 



92$ 2/32 9$ 



9(52 



tl; dp dp 



, 92$ 
9^ 



Ajsm^S 2sin2,5 (25 + 12sin2(5 
1 + sin2 5 



\s\n* 5) 



(1 + sin^ + 4sin^ Sf 







12 



EG. Bizzed, A.M. Bizzed-Sona 



This equation is exactly separable if f3'^Vi{S,/3) does not depend on f3, i.e. 
Vi{6,P) — v{j)/ 0^. This class of potentials has been considered, e.g., by 
Fortunato [24] in his general discussion about the solutions of the Bohr Hamil- 
tonian. We can reasonably suspect that it be not realistic in our case, but a dis- 
cussion of this simply solvable model can shed light on the general properties of 
the quadrupole-octupole vibrations and help to identify more realistic solutions. 
With such a potential, $ — (f){5) is a solution of the differential equation 



A^cj) / , Ajs\i?5 2sin25(25 + 12sin2(5 + 8sin*(5)\ 

■^+\A -«('^) + Y^— 2^ (l+sin25)(l + 4sin2 5)2 j' 



= (26) 



We want to explore first the case of small oscillations of 5 around zero. For our 
discussion, we do not need to specify the exact form of the potential, but, just as 
an example, we can assume a harmonic restoring force and expand up to the sec- 
ond order the terms in sin^ 5 (whose effect, however, will be negligible at least 
for not-to-high values of J), to obtain the eigenvalues A'^ = (2n + 1)Aq and 
eigenfunctions alternatively even or odd, with parity (— )". Or, as an alternative, 
one could use a square-well potential v{6) — Q for 5 < 5t and = +qo elsewhere, 
as exemplified in Fig.|3]a. Other potentials would give a different spectrum of 
eigenvalues but, if they have a single minimum ai 5 — Q, the ground state must 
correspond to a symmetric solution. We indicate with A'_^ the eigenvalue corre- 
sponding to the ground state and with A'_ the one corresponding to the lowest 
antisymmetric solution. We want that, for even spin and parity, the equation in 
(3 take the form of that of the X(5) model, i.e., in our notations. 

Comparison with Eq. |25] shows that this result would be obtained with A'q = 
1/4. Actually, it is probably unnecessary to assume that this relation holds. A 
general problem of all models considering explicitly only part of the overall set 
of dynamical variables, is the effect of the zero-point energies of the neglected 
degrees of freedom, which possibly depend on the value of the active model 
variables. If a phenomenological potential is used for the latter, this potential 
should already include the zero-point energies of all other degrees of freedom 
not explicitly taken into account. 

Assuming that the positive-parity states must be described by the X(5) 
Hamiltonian, the following equation in P holds for the lowest K"^ = 0^ bands 

We now consider in particular the case of the critical-point potential, Vq = 
when (3 is in the interval — /3o and Vq = +oo elsewhere. As for the parity 
dependent term A^, with our assumptions is A-|_ — 0. For negative parity, the 
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Figure 2. The experimental (left) and theoretical energies (normalized to that of the first 
excited state) of the K — bands in ^^"Nd and ^^^Sm. In the latter, the experimental 
values for the K'" — 1~ band are also shown. Theoretical values for even J and parity 
are those of the X(5) model, for odd J and parity are obtained with A_ = 15 for ^^"Nd 
and with A_ = 20 for ^^^Sm. 



term A_ = A'_ — A'_^ will be considered as an adjustable parameter. With this 
assumption, the spectrum of eigenvalues is given by 

e(s, J,^)-eo -C[x(s, J,7r)]2 (29) 

with C constant, x{s, J, tt) the s*'* zero of the Bessel function Jy(a;) and 

ly = v/J(J+ l)/3 + 9/4 + (30) 

Fig Eshows a partial level scheme of ^^^Sm and ^^°Nd, normalized to the en- 
ergy of the first excited state, and compared with the values derived from Eg. 1291 
(which, for positive-parity states, coincide with those of the X(5) model). The 
agreament is fairly good in both cases. In ^^^Sm, the comparison can be ex- 
tended to the lowest negative-parity state of the s = 2 band, if the level reported 
as l'^) really belongs to this band as it would be suggested by its decay. If it is 
so, its energy is significantly lower than the model prediction, but this happens 
also for all the s = 2 excited states of even parity and spin. It must also be noted 
that in ^^^Sm the odd- J levels of the next octupole band {I-C^ = 1^) seem to 
mix with those of the lowest one, as indicated by the large odd-even staggering. 
In the heavier = 90 isotones the second octupole band (with K ^ 0) gets 
closer to the first one, and our model is no longer valid in such a case. In well 
established X(5) nuclei lying in different regions, the lower negative bands seem 
to be built over intruder states, having a band head — 3^ or 5^. 

Is the agreement shown by Fig.|2]a purely accidental one? It is possible, but 
at least it shows that in this case experimental data can be reproduced by using 
a separable potential. We shall see that this is not possible in other cases, and in 
particular in those which seem more interesting to us, ^^^Th and ^^''Ra. 

In these two nuclei, the even J, positive parity states follow closely the X(5) 
model predictions. The negative-parity part of the ground state band lies sig- 
nificantly higher than the positive parity part at low values of J, while the two 
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a 










Figure 3. The lowest even and odd parity wavefunctions in a separable model with 5 
restricted to a small interval around (part a) or around ±5o (part 6). The part c shows 
the corresponding allowed regions in the /32 — /S-j, plane (ci and C2, respectively) and 
the region corresponding to one standard deviation for the localized wavefunctions in the 
schematic model discussed in the text (ca). 



parts merge together above J w 8. This behaviour is very different from those 
of the 7-excited bands in the X(5) nuclei and of the negative parity bands in 
^^^Sm and ^^°Nd. It could be natural to assume that in these mass-224 nuclei 
the axial quadrupole and octupole amplitudes do not oscillate independently, but 
are strictly correlated together': the nuclear deformation oscillates along two 
symmetric valleys, starting at /3 = and caracterized by an average value of S 
close to +So or to —Sq. With reference to our Eq.s|^ this means that /? os- 
cillates in a wide interval, while 6 is confined to a narrow region around some 
nonzero value. As the potential must be symmetric with respect to 6, this means 
that possible values of 5 are localized in one of the two regions around +5^ or 
— ^0- According to Eq.|26l(with indepentent of (3), this would result in a 
negligible parity staggering. 

In fact, as long as we forget the symmetry of the wavefunction with respect 
to (5, we find a doubly degenerate spectrum, corresponding to eigenfunction lo- 
calised either around or around —6q. Reflection symmetry requires that the 
complete wavefunction be symmetric (the sum of the two) for even values of 
J and TT, or antisymmetric (their difference) for odd values. As long as the two 
localized wavefunctions have no overlap (or have a negligible one), the eigenval- 
ues corresponding to even or odd combinations are (almost) equal to those of the 
localized solutions. The dependence on J of the last-but-one term of Eq.l26lhas 
a minor effect on the results, particularly at low values of J. As a consequence, 
the ground-state band has no staggering (or a very limited one). 

How can we proceed? A possible way out is to assume a potential, sym- 
metric in 5, and approaching the harmonic behaviour around 5 k, ±5^. In this 
case, the equation is no longer separable in the entire field. For large values of j3, 
one still finds localized solutions, with a j3 dependent extension around the av- 

' Such a behaviour has been predicted as possible in this region by Nazarewicz and Olanders 
calculations in the frame of the Strutinski model [25]. 
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erage value (±5o): for a fixed value of /3, the variance of 5 would have the form 
D^{d) = C//3. As /3 0, the two localized solutions start to overlap (Fig. |3lc3), 
and this would remove the degeneracy of even and odd combinations, as long 
as the wavefunctions in f3 have a non-negligible value in the region where the 
overlap is significant. At moderately large values of J, the significant values of 
tpif^) are pushed out of this region by the centrifugal-like term of Eq.|28j and the 
staggering will disappear. 

Although this approximation is certainly insufficient just in the region of 
overlap, we think it can give a qualitative explanation of the behaviour of 
the negative-parity band (not far from the point of view of Jolos and von 
Brentano [17]) and, perhaps, give some indication for future developments. 
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